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Abstract 

We present evidence for spin-charge separation in the tunneling spectrum of a system consisting 
of two quantum wires connected by a long narrow tunnel junction at the edge of a GaAs/AlGaAs 
bilayer heterostructure. Multiple excitation velocities are detected in the system by tracing out 
electron spectral peaks in the conductance dependence on the applied voltage, governing the energy 
of tunneled electrons, and the magnetic field, governing the momentum shift along the wires. 
The boundaries of the wires are important and lead to a characteristic interference pattern in 
measurements on short junctions. We show that the experimentally observed modulation of the 
conductance oscillation amplitude as a function of the voltage bias can also be accounted for by 
spin-charge separation of the elementary excitations in the interacting wires. 

Manuscript for proceedings of the NATO workshop on Theory of Quantum Transport in Metallic 
and Hybrid Nanostructures, St. Petersburg, August 2003, to be published by Kluwer Academic 
Publishers B.V., Dordrecht, The Netherlands, edited by V. Kozub and V. Vinokur. Talk was 
presented by B. I. Halperin. 

PACS numbers: 73.21.Hb,71.10.Pm,73.23.Ad,73.50.Jt 



1 



I. INTRODUCTION 



One-dimensional (ID) electronic systems are a very fertile ground for studying the physics 
of interacting many-body systems. Gapless electron gases in one dimension possess universal 
low-energy properties which can be mapped onto the exact solution of the Luttinger model.— 
Such ID systems are collectively termed as Luttinger liquids (LL's). Despite a number of 
remarkable predictions for electronic and thermodynamic properties of LL's,— some made 
more than twenty years ago, direct experimental verification for many of them has remained 
a challenge. This mostly owes to the high quality of quasi-lD systems necessary to bridge 
the gap between the usual Landau Fermi-liquid (FL) physics in three spatial dimensions and 
the LL physics in one dimension. 

An important prediction of LL theory is that the low-energy elementary excitations of 
a one-dimensional metal are not electronic quasiparticles, as in the Landau FL theory of 
three-dimensional Fermi systems, but rather are separate spin and charge excitations that 
propagate at different velocities. An electron entering an LL will split into spin and charge 
excitations, and the electron propagator will have singularities corresponding to both veloc- 
ities, in contrast to the case of a Landau FL where there is a simple pole at a single Fermi 
velocity. In this paper, we discuss evidence for spin-charge separation in tunneling between 
two parallel quantum wires at a cleaved edge of a double-quantum-well heterostructure. We 
use two approaches: one based on mapping out the elementary-excitation dispersions by 
measuring the conductance G as a function of the magnetic field B applied perpendicular 
to the plane connecting the wires and the voltage bias V, and the other focusing on the 
conductance oscillation pattern, in the (V, B) plane, arising due to the finite length of the 
tunnel junction. 

II. EXPERIMENTAL METHOD 

The two parallel ID wires are fabricated by cleaved-edge overgrowth (CEO), see Fig. 
Initially, a GaAs/AlGaAs heterostructure with two closely situated parallel quantum wells 
is grown. The upper quantum well is 20 nm wide, the lower one is 30 nm wide and they 
are separated by a 6 nm AlGaAs barrier about 300 meV high. We use a modulation doping 
sequence that renders only the upper quantum well occupied by a two-dimensional electron 
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gas (2DEG) with a density fiw2x 10 11 cm -2 and mobility /i«3x 10 6 cm 2 V -1 s -1 . After 
cleaving the sample in the molecular beam epitaxy growth chamber and growing a second 
modulation doping sequence on the cleaved edge, two parallel quantum wires are formed in 
the quantum wells along the whole side of the sample. Both wires are tightly confined on 
three sides by atomically smooth planes and on the fourth side by the triangular potential 
formed at the cleaved edge. 

Spanning across the sample are several tungsten top gates of width 2 /iin that lie 2 /j,m 
from each other (two of these are depicted in Fig. [I}. The differential conductance G of the 
wires is measured through indium contacts to the 2DEG straddling the top gates. While 
monitoring G with standard lock-in techniques (we use an excitation of 10 /iV at 14 Hz) at 
T = 0.25 K, we decrease the density of the electrons under a gate by decreasing the voltage 
on it (V^). At V g = V2D, the 2DEG depletes and G drops sharply, because the electrons 
have to scatter into the wires in order to pass under the gate. For V 2 d > V g > Vy the 
conductance drops stepwise each time a mode in the upper wire is depleted.— In this voltage 
range, the contribution of the lower wire to G is negligible because it is separated from 
the upper quantum well by a tunnel barrier. When V g = Vjj, the upper wire depletes and 
only the lower wire can carry electrons under the gate. This last conduction channel finally 
depletes at Vl and G is suppressed to zero. 

The measurements are performed in the configuration depicted in Fig. ^ The source 
is the 2DEG between two gates, g\ and #2 in Fig. ^ the voltages on which are V\ < Vl 
and Vl < V2 < Vu, respectively. The upper wire between these gates is at electrochemical 
equilibrium with the source 2DEG. This side of the circuit is separated by the tunnel junction 
we wish to study from the drain. The drain is the 2DEG to the right of #2 (the semi-infinite 
2DEG in Fig.^J) and it is in equilibrium with the right, semi-infinite, upper wire and with the 
whole semi-infinite lower wire in Fig. Thus, any voltage difference (V") induced between 
the source and the drain drops on the narrow tunnel junction between the gates. In addition, 
we can shift the momentum of the tunneling electrons with a magnetic field B perpendicular 
to the plane defined by the wires. This configuration therefore gives us control over both 
the energy and the momentum of the tunneling electrons. 



3 



III. DISPERSIONS OF ELEMENTARY EXCITATIONS 



The conductance for a spacing of 2 /zm between gates gi and g% is shown in Fig. El The 
measured bright and dark curves in the plot can be interpreted as spectral peaks tracing 
out the dispersions of the elementary excitations in the wires.— In the case of noninteracting 
electrons, the curves are expected to map out parabolas defining the continua of electron- 
hole excitations across the tunnel barrier for various pairs of ID modes, one in the UW and 
the other in the LW. At small voltages, electron repulsion is predicted to split the curves 
into branches with slopes corresponding to different charge- and spin-excitation velocities, 
crossing at a point with V = and magnetic field necessary to compensate for the Fermi 
wave-vector mismatch between the ID modes (in the following referred to as the "crossing 
point" )M 

The ID modes in the upper quantum well are coupled to the 2DEG via elastic scattering 
which ensures an Ohmic contact between the 2D well states and the states confined to the 
cleaved edge. 6 (This scattering is however weak on the scale of the junction length, not 
affecting the finite-size quantum-interference effects in tunneling.™) In addition to tunneling 
between the ID states of the wires, there is a direct electron transfer from the 2DEG to the 
lower wire, when the extended states have an appreciable weight on the edge. Each of the 
quasi-lD wires carries several ID modes. In our theoretical analysis, we will only consider 
the transition between the lowest ID bands of the wires (i.e., the bands with the largest 
Fermi momentum), \ui) «-> which has a strong signal, as seen in Fig. given by a 
family of curves crossing at the lowest magnetic field. Since other ID modes have Fermi 
velocities smaller than \ux) and by at least 40 %,— we can disregard their coupling to 
the lowest bands,— keeping in mind, however, that the effective electron-electron interaction 
is affected by screening due to the nearby 2DEG and other ID modes. Both the spin-orbit 
interaction and Zeeman splitting are negligible in comparison to the Fermi energy, so that 
the electron states are nearly spin degenerate in our heterostructure. In our theoretical 
discussion we therefore consider tunneling between two coupled spinful modes having some 
effective intrawire and interwire interaction. To this end, we use Luttinger-liquid formalism,— 
assuming sufficiently low temperature and voltage bias. It is important in our analysis that 
the measured densities of in the UW and \h) in the LW happen to match to within 
several per cent,— so that the electronic excitations in the double- wire system are collective 
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across the tunnel barrier as well as within each wire. 

The geometry for our theoretical description is shown in Fig. El The potential well U(x) 
is felt by electrons in the upper quantum wire, which are confined to a region of finite length 
by potential gates at both ends (see the source region in Fig. Q). One of these gates (gi) 
causes the electrons in the lower wire to be reflected at one end, but the other (g 2 ) allows 
them to pass freely under it. The effective tunneling region is determined by the length of 
the upper wire, which is approximately the region \x\ < L/2 in Fig. El The magnetic field, 
B, gives a momentum boost hqs = eBd along the x-axis to the electrons tunneling from the 
upper to the lower wire. The current (for electrons with a given spin) 

/OO /■OO 
dxdx' / dte iqB{x ~ x ) e <eVt/R C(x, x'\ t) , (1) 
-oo J — oo 

is determined to lowest order in perturbation theory by the two-point Green function^ 

C(x,x';t) = ([^ u (x,t),^(x',0)]) 

= G>(x, t; x', 0)G<(x', 0; x, t) - G<(x, t; x', 0)G>(x', 0; x, t) . (2) 

The last equality in Eq. (J2J) is valid when the interwire electron-electron interactions vanish. 
Although it might not be a good approximation for our closely-spaced wires, for pedagogical 
reasons we will discuss this limit first. One-particle correlators are defined by the usual 
expressions: G > (x,t;x',t') = -i (v(x, t)^(x', t')) and G<(x, t; x', t' ) = % (&(x', t' )^f(x, t)). 
For V > (V < 0), only the G>Gf {G<Gf) term in Eq. © contributes to the current ©. 

In the LL picture, we can distinguish between the left- and right-moving electronic 
excitations in a given ID mode.— In long wires, each chirality contributes terms propor- 
tional to e ±lk F(x-x') t one-particle Green functions (away from the boundaries), where 
hp = (7r/2)n and n is the electron density in the mode. If the magnetic field is small 
enough, q B kp, the edge-state chirality of the electrons cannot be changed during a 
tunneling event, and the total current is thus a sum of the right-moving and left-moving 
contributions. It is sufficient to calculate the tunneling rate of the right movers only, since it 
equals to that of the left movers under the magnetic-field reversal, B — > —B (so that the to- 
tal conductance is an even function of B). The corresponding zero-temperature one-particle 
Green functions for a gapless translation- and spin-rotation-invariant ID gas of interacting 
electrons has a universal form 

±ip(x)ip* (x 1 ) 



G > ' < (x,t;x',0) 



2w J (z - v a t ± 20+) (z - v c t ± i0+) 



z 2 - (v c t =F ir c ) 2 



(3) 



where v c = Vp/K is the charge-excitation velocity, which is enhanced with respect to the 
Fermi velocity vf by electron repulsion, v s is the spin-excitation velocity, which is close 
to vp for vanishing backscattering rate and is determined by the exchange interaction of 
neighboring electrons for strong repulsion with a sizable backscattering, and a = (K + 
K^ 1 — 2)/4 is a nonuniversal exponent. K < 1 is the compressibility normalized to that 
of the free-electron gas at the same density, z = x — x', and r c is a short-distance cutoff. 
i/j(x) = e ikpX is the noninteracting-electron wave function for an infinite wire at the Fermi 
level. Using Eq. (JSJ) to calculate the two-particle Green function (J2J) and then performing 
the integration in Eq. one shows that the interactions do not shift the position of the 
crossing point at V = and g# = Akp, the mismatch in the Fermi wave vectors of the wires, 
but are manifested by multiple peaks in the tunneling conductance G(V, B), which intersect 
at the crossing point with slopes determined by different spin and charge velocities.— ^ 

It turns out that in a symmetric double-wire structure, the interwire electron-electron 
interactions, do not change the two-point correlation function ((21) apart from renormalizing 
the parameters entering Eq. (jSJ)-^ 1 The reason for this is that a tunneling event from the UW 
to the LW at a low magnetic field creates a long-lived exciton across the tunnel-barrier with 
the electron in the lower and hole in the upper wire, moving in the same direction. The 
exciton propagates freely, as an acoustic plasmon in a single wire but with the velocity v c 
reduced by the electron-hole attraction (corresponding to a larger normalized compressibility 
K). 7 Similarly, v s entering Eq. (J3J) should be thought of as the antisymmetric spin velocity of 
coupled antiferromagnetic Heisenberg chains, again reducing to vf in the case of vanishing 
backscattering. For a symmetric double-wire structure, one therefore expects two velocities 
to be present at low magnetic-field and bias data: the antisymmetric spin- and charge- 
excitation velocities. The situation is more complicated at high magnetic fields capable 
of flipping the electron cleaved-edge chirality: A tunneled particle moves in the direction 
opposite to the hole left behind, decaying into a combination of symmetric and antisymmetric 
excitations across the tunnel barrier, even in the case of a perfectly symmetric double wire.— 
We will not discuss this regime here. 

One can see a family of measured curves crossing at V = and B ~ 0.1 T in Fig. El which 
constitute the signal from the \u\) <-> tunneling. We also draw in Fig. |2]as black solid 
lines the expected parabolic dispersions for noninteracting electrons at the same electron 
densities as \u\) and the white solid lines are generated in a similar way but after 
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rescaling the GaAs band-structure mass, and correspondingly the low-voltage slopes, by a 
factor of 0.7. Remarkably, for positive voltages, we can fit the three visible experimental 
curves, a, b, and c, by such parabolas crossing at V = with two different slopes. Such 
a fitting is a naive extrapolation of LL spectrum at small voltages, where the excitation 
dispersions can be linearized and the two visible low-voltage slopes can be associated with 
elementary ID excitations. Understanding the high-voltage regime, where the dispersions 
acquire a curvature, requires going beyond LL theory, and for that matter beyond the scope 
of this paper. For negative voltages, fewer G(V, B) peaks are visible after we subtract a 
large background signal due to the direct 2DEG-1D tunneling; in particular, only one curve, 
d, is visible for the \ui) «-> tunneling, which cannot be fitted by our naive procedure at 
V < —10 mV. It is important to point out that the observed curves a, b, and c in Fig.|2]rule 
out the noninteracting-electron picture for tunneling between two ID modes with different 
Fermi velocities: If the LW had the higher Fermi velocity, we would expect to see dispersions 
e and a in Fig. |2j and if the UW had the higher velocity, we would see only b and c, but not 
three curves. 

The curve c slope is given by the Fermi velocity of the noninteracting electron gas, 
corresponding to the value of the electron density measured in Ref. [J The slope of a has the 
velocity enhanced by a factor of 1.4. Identifying the faster velocity with the charge mode 
and the slower with the spin mode, we arrive at the following LL parameters characterizing 
the system: K « 0.7 and v s ~ Vp. This value of K indicates that the electron-electron 
interaction energy in the cleaved-edge quantum wires is comparable to the Fermi energy, 
resulting in a sizable effect on the correlation and thermodynamic properties. The closeness 
of the spin velocity v s to the noninteracting Fermi velocity implies a small backscattering 
rate due to electron repulsion, as expected in our wide quantum wires.— 



IV. FINITE-SIZE EFFECTS 



The momentum of the electrons tunneling through a window of finite length L is only 
conserved within an uncertainty of order 27r/L, resulting in conductance oscillations away 
from the main dispersion peaks. We show in this section that another spectacular manifes- 
tation of spin-charge separation at small voltages can be tracked down in such oscillations, 
as those forming checkerboard-like patterns near the crossing points of various dispersion 
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curves in Fig. Ek- We zoom into these oscillations in Fig. 0^. in order to understand them 
in detail, we now generalize our analysis to take into account the finite length of the upper 
wire. 

Assuming that the electron density in each wire varies slowly on the length scale of 
kp 1 (except for unimportant regions very close to the boundaries), we use the WKB wave 
function 

p ikpx p — is(x) 

= 7^= , (4) 

where k(x) = kp[l — U(x)/ Ep] 1 ^ 2 and s(x) = Jq dx'[kp — k(x')}, for right-moving electrons 
in the UW in Eq. U(x) is the potential formed by the top gates defining the finite 
length of the UW, see Fig. El The right movers in the LW are taken to be propagatory, as 
in the infinite wire, assuming its left boundary is formed by the gate gi in the region where 
the UW is already depleted and assuming the gate g 2 potential in the LW is well screened. 
Additional assumptions required for using the WKB wave function (J3J) in the correlator of 
a finite interacting wire are discussed in Ref. Q- 

Substituting Green functions (jHJ) into integral ((TJ), we obtain for the tunneling current at 
V > 

poo 

Joe / dxdx'e i{qB - kF){x - x) i) u {x)i)* u {x')h{x - x') , (5) 



where 

pieVt/h 



/oo 
dt- 
-oo ( 



z 2 - (v c t - ir c y 



(6) 



(z - v s t + iO+)(z - v c t + i0+) 
Here we have taken the two wires to have the same electron density and strength of the 

interactions. For not very strong interactions, a « 1, the last term in Eq. © can be 

disregarded, away from the regime of the zero-bias anomaly— The integrand in Eq. © then 

has two simple poles yielding 

pieVz/(hv s ) pieVz/(hv c ) 

h(z) w -2ttz— — . (7) 

Combining this with Eq. we finally find the differential conductance G = dl/dV: 

G(V, B) oc — — 

where K StC = qs + Akp + eV/(fw s c ) (now including a small mismatch Akp = kp u — kpi in 



-\M(k s )\ 2 - -\M(k c )\ 2 

v., V- 



(8) 



the Fermi wave vectors of the UW and LW, respectively) and 



M{k) = / dx (9) 
<k{x) 



s(x) and k(x) being the same as in Eq. (jlj). M(k) can be found analytically using the 
stationary-phase approximation (SPA): M(k) is evaluated near positions x^ (x + > x~) 
where klx^) = hp — k and the integrand in Eq. Q has a stationary phase. In the case of a 
symmetric potential, U(x) = U(—x) (so that, in particular, x~ = x + ), the SPA gives 

Q(k) r i 

M(k) oc —= v ; cos kx + - s(x + ) -tt/4 , (10) 

JdU(x + )dx 

where Q(k) is the Heaviside step function. The SPA approximation (jlUJl can be shown to 
diverge for small values of k, where we have to resort to a numerical calculation of the 
integral in Eq. JOJ).— The form of Eq. (|10j) shows that (1) the conductance is asymmetric 
in k, vanishing for k < (in the SPA approximation), (2) it oscillates in magnetic field 
(and similarly in voltage) with period Aqs = 2tt/x + , assuming x + to be a slow- varying 
function of k and that (3) it is a superposition of two oscillating patterns in the (V, B) 
plane, the first (second) being constant-valued along line k s = const (k c = const) and 
oscillating perpendicular to it, resulting in a moire structure of G(V,B). [Note that the 
conductance in Eq. (jHJ) so far only includes the right-movers' contribution. In order to get 
the total conductance, one has to add the piece which is mirror symmetric to Eq. (jHJ) around 
B — 0.] If the two velocities v s and v c are not very different, there are two voltage scales 
characterizing the conductance oscillation pattern: 

2nhv c v s , & t r 7ihv c v s , . 

AV= —, -—— r and AKnod = , • H 

AV is the period of the "fast" oscillations, which would be present even in the absence 
of spin-charge separation, and AVmod is the distance between consecutive minima in the 
oscillation power due to the moire amplitude modulation in the voltage direction. The ratio 
between these two scales 

AV mod _ 1 v c + Vg _ 1 1 + v s /v c , . 

AV ~2v c -v s ~21- v s /v c [ } 

can be used to experimentally extract the ratio between the two velocities. We find 

Vs / Vc = 0.67 ±0.07, (13) 

which is independent of the UW length L, while both and AV^d scale roughly as 1/L. 
This value is in agreement with the one found in Sec. IIHI 
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Finally, we compare the interference pattern predicted by our theory, Eq. (JHJ), with the 
experiment, Fig. BJa). G(V,B) calculated using a smooth confining potential given by 
U(x) = Ep exp[(L/2 — |x|)/10] at the boundaries of the upper wire is shown in Fig.^l Many 
pronounced features observed experimentally-the asymmetry of the side lobes, a slow fall- 
off of the oscillation amplitude and period away from the principal peaks, an interference 
modulation along the V-axis, it phase shifts at the oscillation suppression stripes running 
parallel to the field axis-are reproduced by the theory. There is however one experimental 
finding which is not captured by the presented theory: In addition to the periodic modulation 
of the oscillations, there is an appreciable fall-off in amplitude in the voltage direction, as can 
be seen in Fig. EJb). This dephasing can be due to the dispersion curvature which becomes 
appreciable with increasing voltage bias. Its discussion, however, requires going beyond the 
linearized LL theory, which we do not attempt here. 

V. SUMMARY 

The two approaches to study the size of spin-charge separation, one by mapping out the 
dispersions, which are independent of the tunnel-junction length, and the other based on the 
finite-size conductance oscillations, the frequency of which scales linearly with the junction 
length, are found to be in excellent agreement, giving the LL parameter K = vf/v c ~ 0.7 
for the antisymmetric (i.e., excitonic) collective charge excitations in the lowest modes of 
the double-wire structure, and v s ~ vp for the antisymmetric spin velocity. Additional, 
complimentary information about the electron-electron interactions can be extracted by 
measuring the tunneling density-of-states exponent a in the regime of very small voltage 
bias and temperature, where the tunneling rate is suppressed as a power law (the so-called 
LL zero-bias anomaly).^ We do not discuss this in the present paper. 
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FIG. 1: Illustration of the sample and the contacting scheme. The sample is fabricated using the 
CEO method. The parallel ID wires span along the whole cleaved edge (front side in the schematic). 
The upper wire (UW) overlaps the 2DEG, while the lower wire (LW) is separated from them by 
a thin AlGaAs barrier (the wires are shown in dark gray and the 2DEG is light gray). Contacts 
to the wires are made through the 2DEG. Several tungsten top gates can be biased to deplete the 
electrons under them: We show only gi, here biased to deplete the 2DEG and both wires, and 52, 
here biased to deplete only the 2DEG and the upper wire. The magnetic field B is perpendicular 
to the plane defined by the wires. The depicted configuration allows the study of the conductance 
of a tunnel junction between a section of length L of the upper wire and a semi-infinite lower wire. 
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V [mVolt] 



FIG. 2: Measurement of G(V, B) for a 2 /j,ia junction. Light shows positive and dark negative 
differential conductance. A smoothed background has been subtracted to emphasize the spectral 
peaks and the finite-size oscillations. The solid black lines are the expected dispersions of nonin- 
teracting electrons at the same electron densities as the lowest ID bands of the wires, \u\) and 
The white lines are generated in a similar way but after rescaling the GaAs band-structure mass, 
and correspondingly the low-voltage slopes, by a factor of 0.7. Only the lines labeled by a, b, c, 
and d in the plot are found to trace out the visible peaks in G(V, B), with the line d following the 
measured peak only at V > — 10 mV. 
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(b) 



(a) 




FIG. 3: Schematics of the circuit (a) and the model (b). A wire of length L runs parallel to a 
semi-infinite wire. The boundary of the upper wire is formed by potential U(x) confining the one- 
electron wave function ip(x) along the wire. The energy and momentum of the tunneling electrons 
are governed by applied voltage V and magnetic field B. 
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FIG. 4: Nonlinear conductance oscillations at low field from a 6 /im junction, (a) shows the 
oscillations as a function of both B and V. (A smoothed background has been subtracted to 
emphasize the oscillations.) The brightest (and darkest) lines, corresponding to tunneling between 
the lowest modes, break the V-B plain into regions I, II, and III. Additional positively-sloped 
bright and dark lines in II arise from other ID channels in the wires and are disregarded in our 
theoretical analysis. Also present is a slow modulation of the strength of the oscillations along the 
abscissa, (b) Absolute value of the peak of the Fourier transform of the conductance at a fixed V 
in region II as a function of V. Its slow modulation as a function of V is easily discerned. 
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eVL/(27iHv F ) 

FIG. 5: The differential conductance interference pattern near the lower crossing point calculated 
using a smooth confining potential for the upper wire. v c = IAvf, v s = vf, and Afcp = Wir/L. 
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